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ABSTRACT 

Green Bank Telescope (GBT) A21cm observations have revealed a faint, yet 
extensive H I cloud population surrounding the Andromeda Galaxy (M31). The 
newfound objects are hkely analogs to the high-velocity H I clouds (HVCs) seen 
around the Milky Way. At least 20 discrete features are detected within 50 kpc 
of the M31 disk, with radial velocities that are comparable to those of outer disk 
rotation. In addition, a filamentary "halo" component of at least 30 kpc extent 
is concentrated at the M31 systemic velocity. Some of the discrete features are 
organized into elongated systems with velocity continuity, suggestive of tidal 
streams. The discrete population can be characterized by a steep power-law 
distribution of number versus H I mass in the range between 10^ and 10^ Mq. 
The velocity line- width of discrete clouds is correlated with the cloud H I mass: 
such that if the clouds are gravitationally bound this implies a dark- to H I mass 
ratio of ~ 100:1. Possible origins for the discrete and "halo" M31 features include: 
a Local Group "cooling flow", tidal debris from recent mergers or interactions, 
and the gaseous counterparts of low mass dark-matter halos. 

Subject headings: dark matter — galaxies: formation — galaxies: individual 
(M31) 
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1. Introduction 

Much of the sky is covered with neutral hydrogen clouds whose velocities, though within 
several hundreds of km s~^ of zero, are nonetheless anomalous and cannot be explained by 
normal rotation of the Galaxy (MuUer et al. 1963; Wakker & van Woerden 1997). It is 
attractive to consider that these high- velocity HI clouds (HVCs) might be material left over 
from the formation of the Milky Way, now in the process of being accreted (Oort 1970). If 
this is the case, other galaxies may be expected to have similar gaseous components. A key 
advantage of studying such extragalactic clouds is that their distance can be comparatively 
well-constrained. In practice, it has proven difficult to achieve the combination of high 
spatial resolution, mass sensitivity, and large field-of-view necessary to detect HVC analogs 
in external systems. Many nearby galaxies have a few peculiar H I features with H I mass 
as low as ~ 10^ Mq (eg. NGC 628, Kamphuis & Briggs (1992)), but only recently have 
instruments had sufficient sensitivity to begin revealing complete systems of peculiar velocity 
H I (eg. the slowly rotating "halo" component of NGC 2403, Praternah et al. (2002)). In this 
Letter, we report on deep, wide-field observations conducted with the Green Bank Telescope 
(GBT) which have resulted in the detection of a complex system of HI gas extending to at 
least 50 kpc radius around M31. 



2. Observations and Data Analysis 

The GBT observations were obtained over six nights during September 2002. The 
telescope was scanned over a 7° x 7° region around M31 while H I spectra were measured 
over a heliocentric velocity range between —827 km s~^ and -1-226 km at a velocity 
resolution of 1.25 km s~^. At M3rs distance of 770 kpc (Preedman & Madore 1990), the 
GBT beam (9.1' FWHM) projects to 2.0 kpc and our survey covers a 94 x 94 kpc^ region. At 
a resolution of 3 kpc and 18 km s^^ the acheived rms flux sensitivity was 7.2 mJy beam^^ 
corresponding to an H I column density of 2.5 x 10^^ cm~^ in a single channel averaged over 
the beam. 

Our GBT data show many H I clouds surrounding M31 at —515 < Vhei < —172 km s~^, 
which can be compared to the systemic velocity of —300 km s~^. Emission features were 
identified as discrete sources if they were distinct from the M31 H I disk (which extends over 
the velocity interval —620 to —20 km s~^), were unconfused with the Galaxy, and were char- 
acterized by a peak flux density > 5a. Additional GBT spectra were obtained for the faintest 
detections to verify the accuracy of our procedures. We limited the investigation to veloci- 
ties more negative than —160 km s~^ in order to exclude contamination from intermediate 
velocity Galactic emission. 
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3. HI Clouds Around M31 

Our A21cm observations have revealed a population of about twenty discrete clouds 
together with an extended, filamentary "halo" component. The data show that these features 
cannot be instrumental effects related to M31's bright H I disk, for there is no relationship 
between the diffuse structures and the brightness of nearby disk gas. GBT images showing 
the distribution of H I at different velocities are presented in Fig. 1. These were taken from 
a heavily-smoothed variant of the data cube. Contours are drawn at 2, 4, and 6a. Discrete 
compact H I clouds appear at velocities between —172 and —515 km s^^. Extended regions 
of filamentary, low Nhi gas are found preferentially within ~ 80 km s~^ of M31's systemic 
velocity (vhei —300 km s~^). 

We believe that this gas is associated with M31, for the following reasons: (1) there 
is evidence for interaction between several of the clouds and M31 or its companions, (2) 
cloud velocities match the velocity extent of M31's disk and generally correlate with the 
pattern of Andromeda's outer disk rotation, (3) the extended H I component is concentrated 
near M31's systemic velocity, (4) confusion with the foreground Magellanic Stream at this 
location is unlikely, and (5) dynamical constraints on the mass of M31, assuming a bound 
cloud population, are in good agreement with independently determined values. We first 
examine this evidence linking the H I detections to M31, then consider possible origins for 
the cloud population. 

Figure 2 depicts the high velocity gas in the survey field which could be cleanly separated 
from the disk emission of Andromeda. Contours of integrated H I column density at 0.5, 1, 
2, 5, 10 and 20 x 10^^ cm~^ (evaluated using the smoothed cube shown in Fig. 1) are overlaid 
on an optical image of M31. Star symbols mark discrete clouds meeting our 5a significance 
threshold which are adequately resolved from neighbors in position-velocity space. One 
discrete cloud near the north edge of the field peaks below the minimum contour after 
datacube smoothing, but has been independently confirmed in follow-up GBT observations. 
In addition to discrete clouds, there are extended filamentary complexes of low column 
density H I, most notably at small galactocentric radii southeast and north of Andromeda's 
stellar disk. Only structures which are cleanly separated from the M31 disk H I emission are 
included in Fig. 2. Additional filaments can be discerned in Fig. 1. 

The collection of discrete H I clouds extending more than 1.5° (20 kpc) to the south 
of M3rs disk near (aj2ooo; (^J2ooo) = (00:42, 39:15) is suggestive of a tidal stream, based on 
the elongated distribution and continuously varying radial velocity of the features. These 
H I clouds are partially coextensive with the metal-rich stellar "Andromeda Stream" (Ibata 
et al. 2001; Irwin et al. 2002; Ferguson et al. 2002; McConnachie et al. 2003) although they 
are significantly offset to the southwest. The conclusion that these particular clouds are 
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associated with M31 seems inescapable. A second case of partial correspondence of gaseous 
and stellar components is the discrete object which is displaced by only 25' to the southwest 
of NGC 205 and is overlapping in velocity with that galaxy. This feature also merges in 
position and velocity with the M31 disk. 

The only object in our sample that was known previously is Davies' Cloud (Davies 1975), 
which is at least 10 times brighter than the other features. Davies considered that it might 
be part of the Magellanic Stream rather than M31, but in our complimentary wide-field 
study of a 60° x 30° region centered on M31 (Braun et al. (2003); see also Braun & Thilker 
2003, in prep.) we find that the Magellanic Stream comes no closer than 7° to the southwest, 
where its velocity differs by about 100 km s^^ from the diffuse components near M31. It 
seems likely that Davies' Cloud and the other detections from our survey are not confused 
by emission from the Magellanic Stream. A high resolution imaging study of Davies' Cloud 
(de Heij et al. 2002a) has provided additional morphological and kinematic evidence for tidal 
interaction with M31. We suggest that it is simply the most massive of a population of faint 
H I clouds around the Andromeda Galaxy. 

Additional evidence for association of the newly-detected clouds with M31 is that the 
discrete cloud velocities appear partially correlated with the pattern of outer disk rotation in 
M31, with the most negative velocities occurring in the south-west and the most positive in 
the north-east. The high negative velocities of some clouds (e.g., —515 km s~^) are also more 
extreme than any seen previously in Galactic HVC surveys (v > —466 km s~^) (Putman 
et al. 2002; de Heij et al. 2002b). A straightforward interpretation is that the clouds are 
collectively under the gravitational influence of the Andromeda Galaxy and have interacted 
dissipatively with M31 in the past. 

If the detected clouds represent a circum-galactic population gravitationally bound to 
Andromeda, their velocity offsets and projected separations should be collectively related 
to M3rs total mass. Using the virial theorem and representing the discrete GBT clouds 
as bound test particles orbiting a central mass (Binney & Tremaine 1987, Eqn 10-22), we 
find that M^si = 5.2 x 10^^ Mq. This estimate reflects only the mass contained within the 
radii probed by our clouds. If we instead assume an extended mass distribution (Eqn 10-23) 
similar to the distribution of clouds, the entire virial mass of Andromeda is then MM3i,vir 
= 7.9 X 10^^ Mq. These two mass determinations are in good agreement with independent 
total mass estimates of 4 x 10^^ Mq within a radius of 30 kpc (Brinks & Burton 1984) and 
1.6 X 10^^ M0 within Andromeda's virial radius of 300 kpc (Klypin et al. 2002). 

From all these considerations, we consider it very likely that the newfound H I clouds 
reside in the circum-galactic environment of M31. 
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4. Discussion 

At a distance of 770 kpc, most clouds in our sample have H I masses in the range 
0.15 — 1.3 X 10^ M0. Such clouds would not have been detected as discrete objects in 
previous surveys of more distant galaxies, which generally have probed only the mass range 
above ~ 10^ Mq (Braun & Burton 2001). An indication for the more general existence 
of low mass circum-galactic populations comes from the recent detection of a 15% excess 
of H I mass found on scales of several hundred kpc relative to that seen at tens of kpc in 
the environment of apparently isolated galaxies (Braun et al. 2003). Placing accurate limits 
on the total HI mass of the M31 halo cloud population is difficult because we lack a well- 
measured radial scale length for the population, must cope with potential incompleteness 
due to Galactic confusion, and have no conffimed minimum cloud mass. To the degree that 
the observed cloud mass function is consistent with a sensitivity-bounded power law given 
by dN{Miii) a M^f dMni, our estimate of HI mass for the halo cloud population within 
the GBT field is ~ 3 — 4 x 10^ Mq. The upper limit of this range is calculated assuming 
the true population is bounded by a minimum HI cloud mass of 10^ Mq and that there are 
fewer than ~ 100 of such objects, the vast majority of which are not yet detected. We also 
included an extra 25% to account for clouds confused with Galactic emission, but even so, 
the H I mass of Andromeda's discrete halo cloud population in our data amounts to only 1% 
of the mass of M31's HI disk. Below, we explore the idea that the halo clouds may trace a 
more substantial amount of ionized gas and possibly also dark matter. 

An obvious source of high velocity H I is tidal stripping from recent or ongoing mergers. 
We detect a gaseous feature which is partially coextensive with the stellar Andromeda Stream 
(McConnachie et al. 2003), and of comparable spatial extent, as well as a component of 
possibly stripped gas adjacent to NGC 205. However, most of the discrete clouds detected in 
our GBT survey are rather isolated, lacking any apparent relation to known M31 companion 
galaxies, and do not have obvious indications if internal tidal distortions. Nevertheless, the 
faintest dwarf galaxies and stellar streams are notoriously hard to detect (Armandroff et al. 
1999). One might argue that some dwarfs have yet to be discovered, while others have been 
entirely consumed by M31. Indeed, the comparatively young (6-8 Gyr) halo stars found by 
Brown et al. (2003) in the ACS Andromeda deep field do suggest that M31 has seen a major 
merger (with a massive star forming galaxy) or several minor mergers. 

Goohng of a tenuous inter-galactic medium (Oort 1970) is a second viable source of 
clouds. In this view, halo clouds would condense from, and remain confined by, coronal 
gas which is located around M31 or perhaps pervades the Local Group. Sembach et al. 
(2002) presented evidence for an extended and highly ionized Galactic corona or Local Group 
medium traced by high-velocity O VI absorption. Likewise, Lockman et al. (2002) conducted 
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a sensitive HI survey of 860 sightlines at 5 > —43°, which suggested that a low-column 
density "mist" of high-velocity neutral gas surrounds the Milky Way (Nhi ^ 8 x 10^^ cm~^ 
along 40% of sightlines). Together, these studies show that the Galactic HVC phenomenon 
extends to much lower column densities than traditionally appreciated and that the classical 
Galactic HVCs are only part of a more ubiquitous multi-phase medium. The filamentary 
"H I halo" we detect concentrated on the M31 systemic velocity may be a manifestation of 
an M31 "coohng flow" . In such a scenario, the detected H I may represent only the tip of 
the iceberg, in terms of baryonic mass, since the gas may well be only of order 1% neutral 
at the relevant low volume densities (Sembach et al. 1999). The primary source of such 
coronal gas may actually be the action of a large-scale "galactic fountain" (Shapiro & Field 
1976; Bregman 1980; Savage et al. 2003) in the recent or distant past. 

A third component of high velocity H I near M31 might be the gas associated with a 
putative population of low mass dark-matter halos. Current simulations of the Local Group 
in a ACDM cosmology predict a large population of low mass DM halos at the present epoch 
(Klypin et al. 1999; Moore et al. 1999) which dramatically outnumber known dwarf galaxies. 
If the M31 halo clouds are tracers of substantial dark-matter concentrations, this should 
be reflected in their internal line-widths. We plot the observed H I mass versus FWHM 
line-width for discrete clouds in Figure 3. Line-width measurements were determined using 
Gaussian fits to cloud spectra. The line- width distribution has an approximate lower bound 
roughly consistent with thermal broadening for gas at 10^ K (FWHM = 24 km s~^). A 
systematic increase in FWHM line-width is observed with increasing H I mass, reaching 
some 70 km s"^ for the clouds near 10^ Mq. Although Davies' Cloud is significantly offset 
from the remainder of the distribution in Fig. 3, high-resolution imaging has revealed much 
higher internal line- widths in that object as well (de Heij et al. 2002a). To demonstrate 
the expected distribution of line-width with mass, in Fig. 3 we show a curve corresponding 
to = 100GMhiR~^, where the characteristic discrete cloud radius, R, has been held 
constant at 500 pc, based on our interferometric WSRT detections of some of the cloud 
cores (Braun et al. 2003, in prep.). This curve, corresponding to a dark to H I mass ratio of 
about 100:1, is not intended to fit to the data, but merely to provide a basis for comparison. 
The hypothesis of a kinematically dominant dark- matter component appears to be consistent 
with the observed hne-widths of the discrete M31 halo clouds. 

Another aspect of the dark- matter mini-halo scenario that can be checked is the expected 
number of such objects in the appropriate mass range in the vicinity of M31. Sternberg et 
al. (2002) predict ~ 25 dark matter mini-halos associated with gravitationally confined HI 
within a radius of 40 kpc around M31, based on the simulations of Klypin et al. (1999) 
and Moore et al. (1999). This is fully consistent with the 20 discrete M31 halo clouds we 
have detected. The Sternberg et al. calculations also suggest that circum-galactic objects 
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of such low mass and peak H I column density should be only ~ 10 % neutral, implying 

Mhi+hii ~ IOMhi and Mjjm ~ ^OMhi+hii- Associated ionized gas could perhaps be 
detected via deep Fabry-Perot imaging (Bland-Hawthorn & Maloney 1999; Tufte et al. 2002) 
or absorption line spectroscopy of background quasars (Murphy et al. 2000; Tripp et al. 2003; 
Sembach et al. 2003). 

In summary, our GBT observations of the Andromeda Galaxy have revealed the first 
extensive extragalactic counterpart to the Galactic HVC population. Both discrete and 
diffuse components are detected. We find some supporting evidence for at least three different 
possible origins of the high velocity gas, namely: tidal disruption, halo condensation and 
association with low-mass DM halos. 
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Fig. 1. — The overall distribution of high- velocity H I gas in the GBT survey volume, 
displayed with a square-root transfer function (0 — 4.5 Jy beam~^). Contours are drawn at 
2, 4, and 6a. Our data have been smoothed to (3 kpc, 72 km s~^) resolution for this figure. 
Note the population of discrete high-velocity clouds, plus an extended distribution of gas 
present in the channels near M31's systemic velocity (-300 km s~^). Images for the following 
hehocentric velocities are shown: -656 (top left), -622 ... -181 km s~^(bottom right). 
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Fig. 2. — Total column density for discrete and diffuse high- velocity H 1 in the M31 GBT field, 
after masking emission from Andromeda's inclined, rotating disk. Contours were evaluated at 
(3 kpc, 72 km s~^) resolution and rendered at 0.5, 1, 2, 10, and 20 x 10^^ cm~^, then overlaid 
on a V band image of M31. The position of each discrete cloud (or distinct substructure) 
is marked using a star symbol. We indicate the location of the metal-rich stellar stream 
analyzed by Ferguson et al. (2002), by plotting the imaging fields of McConnachie et al. 
(2003). H I emission from NGC205 has been masked. 
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Fig. 3. — Observed H I mass versus FWHM line-width for discrete clouds near M31. Line- 
width was measured using Gaussian fits to flux-weighted average spectra determined over 
the extent of each cloud. We have also plotted the anticipated relationship between H I mass 
and line-width for gravitationally confined clouds. To first order, our FWHM measurements 
appear consistent with the hypothesis that the objects arc dark matter dominated, assuming 
a dark matter to H I mass ratio of 100:1 and a characteristic size of 0.5 kpc for each H I 
cloud core. The later assumption is supported by our WSRT detection of resolved, high 
column-density cores {Nhi ~ 10^^ — lO^'^ cm~^) within many of the more centrally-located 
clouds. 



